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Abstract 

Non-relativistic potential models are considered of the pure power V(r) = sgn(q)r q 
and logarithmic V(r) = ln(r) types. It is shown that, from the spectral viewpoint, 
these potentials are actually in a single family. The log spectra can be obtained from 
the power spectra by the limit q — > taken in a smooth representation P n e(q) for 
the eigenvalues E n t{q). A simple approximation formula is developed which yields the 
first thirty eigenvalues with error < 0.04%. 

PACS 12.39.Pn: 12.39.Jh; 03.65. Ge. 
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1. Introduction 

We consider a single particle that moves in a central potential V(r) and obeys 
non-relativistic quantum mechanics. We study two cases: (a) V(r) = sgn(q)r q , and (b) 
V(r) = ln(r). The main purpose of the paper is to show that these two problems are 
intimately related, and to provide a single formula which yields accurate approxi- 
mations for the corresponding discrete Schrodinger eigenvalues. The power-law and 
logarithmic potentials continue to be employed as non-relativistic models for quark 
confinement [1-5]. Hence it is important to understand that they are, in a sense, from 
the same family of potentials. A detailed analysis of the properties of non-relativistic 
potential models is given in the review article by Quigg and Rosner [6] . The approach 
to this question based on the expression {dr g /dq} q ^ — lnr is also been discussed in 
Ref. [7]. In the present paper we approach the same problem by using what we call 
the ' P -representation' for the discrete spectrum generated by power-law potentials. 
This allows us to construct a general energy eigenvalue formula which joins q = -1 
smoothly to q = 2, passing through the logarithmic case q = 0. 

The power-law family and logarithmic potential will be kept distinct until it 
becomes clear in what sense the logarithmic potential corresponds to the limit q — > 
of the power-law class. We first treat the elementary issue of scaling. We write the 
eigenvalues of the bare problems as: 



where n = 1,2,3, ... counts the discrete eigenvalues in each angular-momentum sub- 
space labelled by £ = 0, 1, 2, . . . . The eigenvalues so labelled have degeneracy precisely 
21 + 1. Scaling arguments [8,3] can be used to show that more general Hamiltonians 
have the corresponding eigenvlues given in terms of {E nt {q), E^ e } by 



where v > is a coupling parameter and ^ > might, for example, be fi = h 2 /(2m). 
As a consequence of these scaling rules we shall only need to consider the special case 
li = v = 1 for the rest of the paper. 

In Section (2) we introduce an exact semi-classical representation P n e(q) for these 
eigenvalues which will allow us to show how the log spectrum can be derived as a 




(1.1) 




(1.2) 
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limit from the family of power spectra. As a motivation for this, we plot the power 
eigenvalues E ne (q) as functions of q in Fig.(l). These graphs are to be compared 
with the corresponding smooth P -representation shown in Fig. (2). Once the P- 
representation is established, we go on to consider the log-power theorem P n t(0) = P£ e , 
and, in Section (3), to develop an approximation formula for P n e(q)- 

2. The functions E(q) and P{q) 

The P -representation for the spectrum is as follows: 



Pi 

r>0 I r 2 



E ni =mm\^f + V(r)\. (2.1) 



This transformation is well defined provided we can always find the P which cor- 
responds to a given E. For smooth potentials, the minimum exists and is unique if 
r 3 V'(r) is monotone, a condition which is certainly satisfied if V{r) is a pure power or 
log. For these potentials, Eq.(2.1) is therefore an exact semi-classical representation 
for the eigenvalues: the kinetic energy term scales like L~ 2 , as it should. It is perhaps 
interesting also to note that positive factors such as n and v can be re-introduced 
in front of the kinetic- and potential-energy terms, without having to revise the P 
numbers. For the specific pure-power and log problems at hand, we have respectively: 

(2.2) 



E ni (q) = mm <{ + sgn(q)r q 

- g »(,)(i +1 )(f 



0p2 \ q+2 



and 

r>0 



^=min^^+ln(r) 



(2.3) 



= ln((2e)ip^). 

Some of the E, and the corresponding P, are already known exactly, from elemen- 
tary quantum mechanics. Thus from the known eigenvalues for the Hydrogen atom 
E n t(—1) = -[2(n + l)]~ 2 and the harmonic oscillator E ne (2) = An + 2£- 1 we immediately 
obtain the outer P numbers of our range of interest -1 < q < 2 : 

P nt {-l)=n + t, P ne (2) = 2n + £ - ^. (2.4) 

Our goal is now to interpolate P n e(q) at interior points of the range. 

The spectral relationship of pairs of potentials {U, V} that are smooth transfor- 
mations V(r) = g(U(r)) of one another may be analysed according to the following geo- 
metrical reasoning [9]. We suppose that the eigenvalues of = -A + vll{r) are given 
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by E [ n /(v). Each tangent to V = g(U) is a potential of the form V*(r) = a(t)U(r) +(3(t), 
where t is the point of contact, and a(t) = g'(U(t)) and /3(t) = g{U(t)) - U(t)g'(U(t)). 
When g has definite convexity, every tangential potential lies on one side of V and 
we obtain energy bounds; an optimization over t yields the best such bound. In the 
concave case ( g" < 0), for example, we obtain: 

EZ\v) < rmn{E^(va(t))+(3(t)} • (2-5) 

The analytical complexity of the optimization in (2.5) is greatly reduced by the in- 
troduction of 'kinetic potentials' [9], which, in the case of power- law potentials, lead, 
after a change of variable and the correct choice of P, to the very simple form (2.1). 
It is straightforward to show for the present problem that each power potential with 
power qi is a convex transformation of all other power potentials with powers q 2 < qi. 
Meanwhile, the log potential is a convex function of each negative power, and, at the 
same time, it is a concave transformation of each positive power. Two results which 
follow from this analysis are (i) [8] the functions P n i{q) are monotone increasing in q, 
and (ii) [9] the log-power theorem, which says that lim g ^ P n e(q) = P n e(0) = P„e- 

As is clear from Fig.(l), the spectra E ni (q) of the power potentials converge on 
the values ±1 as the power q approaches zero respectively from positive or negative 
values. Meanwhile, the set of points for the log spectrum E% e would appear to be quite 
unrelated, and would have no natural place on this graph. In the P -representation the 
picture is quite different, as is shown in Fig. (2). The monotone P(q) curves are distinct 
at q = 0, and we know that they are exactly equal there to the values corresponding 
to the log potential. This smooth expression of the entire spectral family suggests 
that we can interpolate the P curves with a simple polynomial, a task which is taken 
up in the next Section. We made an earlier attempt at an interpolation [8], which was 
much less accurate because it could not benefit from the (yet undiscovered) log-power 
theorem. 

3. Interpolation formula for P(q) 

As a model for the P curves we use a cubic polynomial which allows us to fit values 
at the four points q = {-1,0,1,2}; the exact formulas (2.4) are used for q = {-1,2}. 
We make the expansion about q = so as to favour the q values between log and 
linear that are used in potential models. Thus we define 



P(q) =a + bq + cq 2 +dq 3 , 



(3.1) 
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in which, for simplicity, we have omitted the quantum-number subscripts nl. By 
inversion we have 
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(3.2). 



For the linear potential g = 1 the exact S-state eigenvalues may be expressed in terms 
of the zeros of the Airy function [10]. We take these known values and complement 
them with others computed numerically for the log and linear potentials. By inverting 
Eq.(2.2) and Eq.(2.3) we have thus computed the coefficients P(0) = log and P(l). 
The P-data for the first 30 eigenvalues are exhibited in Table(l), along with the 
approximations we get for the energies E ni {\) For comparison, we have tabulated 
also the percentage errors in these eigenvalue results. For the first 30 eigenvalues, 
the percentage errors are all positive and less than 0.04% at q = \-. as q approaches 
q = or q = 1, these errors decrease dramatically. Similar accuracy is obtained for 
the whole range -1 < q < 2 of the interpolation. 



4. Conclusion 

In spite of present-day computing convenience, it is still very useful to have an 
approximate eigenvalue formula, particularly a simple and accurate one. The possi- 
bility of a simple formula is a consequence of the existence of the smooth monotone 
P -representation for the eigenvalues. Since the early work on potential models there 
has been a realization that, from a practical point of view, a log and a power po- 
tential, with < q < 1, could be adjusted to serve in the model almost equally well, 
especially if q is small. The log-power theorem provides a theoretical basis for these 
concrete spectral observations. 
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Table 1 The 'input' values of P n e(0) and P n e(l); and the approximations E^ e {\) for 
E n t{\) obtained via the cubic P formula (3.1), with the percentage errors. 
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Figure 1. 

The first 30 eigenvalues E nt {q), 1 < n < 5, < ^ < 5, corresponding to the power 
potential V(r) = sgn(q)r q . For q > 0, the eigenvalues increase with q from 1 to 
E nt {2) = 4n+2f+l; for q < 0, they decrease (as q increases) from E ne (-1) = -{2(n+l)}- 2 . 
to -1. Both sets of curves increase with n and I. 
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Figure 2. 

In the P -representation, the same set of 30 eigenvalues shown in Fig.(l) now lie on 
monotone smooth curves. The log-power theorem states that the P values for the 
log potential are precisely P n e(0)- As q increases from -1 to 2, the degeneracy of 
the Coulomb problem P n e(-1) = n + £ evolves into the degeneracy of the harmonic 
oscillator P ne (2) = 2n + P.-\. 



